Quasi ID Bose-Einstein condensate flow past a nonlinear barrier 
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The problem of a quasi ID repulsive BEC flow past through a nonlinear barrier is investigated. 
Two types of nonlinear barriers are considered, wide and short range ones. Steady state solutions 
for the BEC moving through a wide repulsive barrier and critical velocities have been found using 
hydrodynamical approach to the ID Gross-Pitaevskii equation. It is shown that in contrast to 
the linear barrier case, for a wide nonlinear barrier an interval of velocities < v < v~ always 
exists, where the flow is superfluid regardless of the barrier potential strength. For the case of the 8 
function-like barrier, below a critical velocity two steady solutions exist, stable and unstable one. An 
unstable solution is shown to decay into a gray soliton moving upstream and a stable solution. The 
decay is accompanied by a dispersive shock wave propagating downstream in front of the barrier. 

PACS numbers: 42.65.-k, 42.50.Ar, 42. 81. Dp 
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I. INTRODUCTION 

The problem of the transcritical flow of a BEC through 
the penetrable barriers has been under recent active in- 
vestigations [H-i, HJ. The damping processes for the 
superfluid flow moving through the barrier are of a fun- 
damental interest. In multidimensional case above some 
critical velocity of the obstacle motion the damping ac- 
companied by the radiation emission [2] is observed. 
Thus in the region when the motion is still superfluid, 
the velocity is bounded above. The damping is associ- 
ated with the Landau type damping and related to the 
emission of the elementary excitations. Landau damping 
can be described in the framework of the mean field the- 
ory and is not associated with thermalization processes 
Q. The critical velocity value at which the damping is 
observed, differs essentially from the values predicted by 
the Landau theory. As it was shown firstly by Feynman 
, the reason is in the nonlinearity of the system. In the 
case of a quasi ID Bosc-Einstein condensate flow, when 
passing through a penetrable barrier, some interval of ve- 
locities w_ < v < exists, where trains of dark solitons 
are generated, that leads to deviation from predictions 
based on the matching with the spectrum of elementary 
linear excitations [l|,|8[].In addition in this range of veloci- 
ties, generation of dispersive shock waves occurs. Exper- 
imental proof of the existence of the velocities interval 
was given in the work 0] . Hakim [§] has indicated that 
for supersonic velocities (including ones above supercrit- 
ical velocity v+) some radiation is still nonzero and its 
amplitude rapidly decreases at the ratio of the potential 
variation length to the GPE coherence length. The am- 
plitude of the wake can be characterized by the Fourier 
transform of the obstacle potential Thus, wide and 
smooth potentials can be considered as radiationless at 
velocities above supercritical. Seemingly in one dimen- 



sional case only stable dark solitons can exist. Peculiar- 
ity of the one dimension is in the fact that generation 
of the solitons is possible till some supercritical velocity, 
v + . Above this velocity the emission is strongly damped 
and the quasi-superfuidity is restored. The radiation ex- 
ists, but exponentially small-decay rate is proportional 
to Ipot/lfn where lh is the healing length of the order of 
the dark soliton width. 

In this work we consider the phenomena occurring in 
the flow of a quasi ID BEC past a nonlinear barrier which 
is a localized space inhomogeneity of the the nonlinearity 
coefficient in the Gross-Pitaevskii equation. Such a type 
of barriers can be formed by some area of BEC where the 
effective value of the atomic scattering length is varied in 
the space. It can be achieved both by the Feshbach res- 
onance techniques [ll|, and by the local variation of the 
transverse frequency of the trap potential. In the former 
case, varying external magnetic field in space near the 
resonance, one can vary the value of the atomic scatter- 
ing length a s . Another way is to use optically induced 
Feshbach resonances [l2j. In this case the variation can 
be achieved by local change in the intensity of a laser 
field. Variation of a s in a half space recently has been 
suggested to generate vortices in BEC as a nonlinear pis- 
ton method [14| . 

The present paper is motivated by the works [l|, Q 
where flow of a BEC past an obstacle in one dimension 
was investigated. We consider two cases, wide obstacle 
potential and short range one. 



II. THE MODEL 

Let us consider a nonlinear penetrable barrier mov- 
ing through the elongated BEC. A quasi one dimen- 
sional BEC can be described by the Gross-Pitaevsky 
(GP) equation with standard dimensionless variables 
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where 



t = Tuj x , x = X/l±, 4>(x,t) = v^kl*(a:,*). 

l x = y/K/muZ, (2) 

a s is the atomic scattering length, uj± is the transverse 
frequency of the trap, V — > a s o is the background 
value of the scattering length a s . For the further study 
of the flow problem it is useful to pass to the reference 
frame moving with the barrier x 1 = x + vt, t = t. So we 
come to the equation 



uu x + p x + (Vp) x = 0, 
with the boundary conditions 

p 1, u v, V(x) — > 0, when \x\ 
Integrating over x we find 

pu = v, 



-u 2 + p + V(x)p = t;v 2 + 1. 



CO. 



(8) 
(9) 
(10) 
(11) 



x' x' 



l^l 2 ^ = V(x')\ip\ 2 i[). (3) Eliminating the function p from these equations, we get 



The scattering length can be manipulated with a laser 
field tuned near a photo association transition, e.g., close 
to the resonance of one of the bound p levels of the ex- 
cited molecules. Virtual radiative transitions of a pair of 
interacting atoms to this level can change the value and 
even reverse the sign of the scattering length [12]. Re- 
cently spatial modulations of the atomic scattering length 
by the optical Feshbach resonance method was realized 
experimentally in BEC [l3| . Such approach implies some 
spontaneous emission loss which is inherent in the opti- 
cal Feshbach resonance technique. Here we assume that 
such dissipative effects can be ignored, since they be- 
come possible if one uses laser fields of sufficiently high 
intensity detuned from the resonance. Thus the repulsive 
nonlinear barrier can be formed by an focused external 
laser beam with the parameters lying near the optically 
induced Feshbach resonance. 



A. Wide obstacle potential 

We analyze this case following the method developed 
in [l], [f| for the linear barrier case. Let us pass to the 
hydrodynamical form for the GP equation (TTJ) . It can be 
obtained by the following transformation 



ip(x',t) = y/ P {x',t)e 1 ^' u(x ' t)dx . 

Substituting it into jl]) and introducing u' 
obtain the system 



(4) 



u + v we 



Pt + (pu') X ' = 0, 

,2 



Px 



4p 



/V + {V(x')p) x , = 0. 



(5) 



(6) 



For a wide smooth obstacle potential we can neglect the 
terms in the bracket in the second equation that corre- 
sponds to the hydrodynamical approximation. Omitting 
also primes, for stationary solutions we can put pt = 
and u t = 0, and obtain the following system of equations 



V(x) = —(u-v)[2-u(u + v)] = F(u). (12) 
2v 

Since we consider repulsive obstacle potential V(x) > 
we have the condition F(u) > 0. Maximum of F(u) is 
realized at u m = \J (v 2 + 2)/3. Thus the maximum of 
the function F(u) is 



max[_F(u)] = /i(v) 




1. 



(13) 



Stationary solution u(x) is obtained by solving the 
equation (TT2"j) with respect to u. This equation has a 
real solution defined for all x provided that 



V m = ma,x[V(x)] < max[F(u)], 



(14) 



i.e. the range of values of V(x), which is [0,1/m], lies 
within the range of values of the function F(u) [l|. 
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FIG. 1: Maximum of the function F(u) (see Eq.([T3]l) 

versus 

x. For given obstacle potential maximum V m = 0.5, critical 
values of the velocity v- — 0.409, «+ = 2.117 



Maximum of the function F(u) versus the obstacle 
velocity v of BEC is presented in Fig. [1] As seen for 
any value of V m two critical values of the velocity exist, 
V—, v + , determined by equation V m = p{v). In transcrit- 
ical regime, in the interval t>_ < v < v +1 the condition 
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of the stationary flow ([14]) does not hold. Out of this 
region, in sub-critical (v < w_) and supercritical (v > v + ) 
regimes the radiation phenomena are negligible and the 
motion of the system can be considered as superfluid. 

Analyzing expression (fl"3f and Fig.[T]it should be noted 
that unlike the case of a wide linear barrier, considered 
in [l| , the velocity u_ is not vanish and there always exists 
an interval < v < u_ where the flow is superfluid. 

Eq. (fT2]) can be rewritten as 



(v 2 + 2)u + 2v(V(x) + 1) = 0, 



(15) 



which is a cubic equation with respect to u(x). Solving it 
we obtain the following solutions for u(x) satisfying the 
boundary conditions 



u(x) = — 2 v / gcos I s(x) — 



u(x) = —2^008 ( s(x) + — 



2n 



2tt 



for v < v~, (16) 
for v + < v, (17) 



where 



q = 



— i s ( x ) = 3 arccos( 



v(V(x) + 1) 



)• 



Spatial profiles of the local velocity u for sub-critical 
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FIG. 2: Spatial profiles of the local velocity u(x). The barrier 
velocities are equal to v = 0.373 and v — 2.517 for lower and 
upper lines respectively. 

v = 0.373 (v < u_) and supercritical v — 2.517 (y > w+) 
regimes are depicted in Fig. [5] The NL obstacle poten- 
tial is taken in the form V(x) = V m / cosh(a;/2) with its 
maximum value V m = 0.5. 

Fig. |3] depicts time evolution of a BEC flow through 
a repulsive non-linear potential V(x) = V m / cosh(x/2) 
with V m = 0.5 in (a) subcritical (v = 0.373 < u_) and 
(b) supercritical (v — 2.517 > v + ) regimes, respectively. 
Initial form of the condensate density p(x) is determined 
by Eq. (|10[) as p(x) = v/u(x), where initial distribution 
of local velocities u(x) is given by Eqs. (|TT))) . (JTTJ) . One 
can see that in these regimes the flow through the barrier 



FIG. 3: Time evolution of a BEC flow in (a) subcritical 
regime, v = 0.373 and (b) supercritical regime, v — 2.517 
through a nonlinear repulsive potential barrier V m / cosh(x/2) 
with V m — 0.5. The initial wave packet and distribution of the 
BEC local velocities u(x) are taken in the form determined 
by formulas (TI5J, (TTJl and Eq. (flu]) . 



is steady. Existence of small amplitude waves, spreading 
from the hump in the beginning is a result of neglecting 
small terms in the course of derivation of Eqs. ([7]) and 
©. In Fig. [5b one can see that in supercritical regime the 
solution at the center has the hump form. The numerical 
simulations show stability of this kind of steady flows. 

In order to carry out numerical simulations of the be- 
havior of a BEC at transcritical velocities (v_ < v < v+), 
we can not use Eqs. (|16p . (JTTJ) as an initial wave packets, 
because they have been derived for a steady flow. 

In numerical simulations it is more convenient to in- 
crease adiabatically the strength of NL potential V m . 
In Fig. [4] we show time evolution of BEC flow through 
a NL potential barrier in the transcritical regime with 
v = 0.47 (v > i>_). The NL potential is taken in the form 
V(x) — V m / cosh(ir/2). V m is increasing from to 0.5 in 
the time interval < t < 1000 and then is kept constant. 
One can see that in the transcritical regime the flow be- 
comes unsteady and a train of dark solitons emerges from 
the NL barrier at the barrier potential strength V m = 0.5. 



4 




110 



where unknown parameter xq depending on the poten- 
tial strength 7 is determined from the relation 



7 = 



(1 — v 2 ) 3 / 2 cosh(Vl — v 2 xo) sinh(\/l — v 2 xq) 
(v 2 + sinh 2 (\/l - v 2 x )) 2 



(23) 



obtained from matching condition for derivatives R x (x) 
at x = 

R x (+0) - R x (-0) =7^ 3 (0). 



FIG. 4: Time evolution of a BEC flow in the transcritical 
regime when the NL barrier velocity v — 0.47 (v- < v < 
«+). The NL barrier is taken in the form of V m / cosh(a;/2) 
with Vm = 0.5. During the time period from t = to t = 
1 000 (that is not presented in the figure) the value of V m is 
adiabatically being increased from to 0.5. Further evolution 
is given at V m = 0.5. 



B. Short range nonlinear obstacle (delta- function 
potential) 

In this section we follow the approach used in the work 
Q . Let us suppose the condensate to have a chemical po- 
tential ji = 1. Then in the frame of the moving obstacle 
with the velocity v equation ((T|) takes the form 

ii/H + iv^x + ^» - V - = V(x)\i/j\ 2 ip (18) 

with uniform boundary conditions ^(a;)! 2 = 1 at i -> 
±00). 

Looking for time independent solution in the form 
ij)(x) = R(x) exp(i(f>(x)) we get equations for amplitude 
R(x) and phase cf>(x) 



> • I 1 - jp ) > ( 19 ) 



Rxx = v 2 [-R + ^3 J + R 6 + V{x)R i - R. (20) 

In the case of the 5 function barrier potential (a sharp 
jump in the nonlinearity) V(x) = j5(x) the solution R(x) 
has the form 



R 2 (x) = 1 



1 - v 2 



cosh [VT — v 2 (x =F xo)] 



at x § 0, (21) 



Substituting obtained R(x) into Eq. (fT9|) and solving it 
we obtain phase <fi(x) as 



arctan 



2wVl - v 2 



4>(x) = f{x) = 

at x > 0, 



exp(Vl — v 2 (x + xq)) + 2v 2 — 1 

and 4>(x) = 2/(0) - f(-x) at x < 0, (22) 




FIG. 5: Dependence of the parameter xq on the nonlinear 
potential strength 7 for f = 0.65. 

Fig. [5] depicts a typical relation between the potential 
strength 7 and parameter xo at v = 0.65. As seen for 
given strength 7 there are two values of the parameter 
xq (or not a single) corresponding to a pair of steady 
solutions. One of the solutions (xn = xq—) is unstable 
and another (xo = x +) is stable 0, @. 

Time evolution of stable and unstable steady solutions 
corresponding to x$ + = 0.752048 and xq- = 0.350966 
are shown in Fig. [6] As seen the unstable solution decays 
into a gray soliton moving upstream with the velocity less 
than v and a stable solution localized at the barrier posi- 
tion. The decay is accompanied by the radiation emitted 
downstream in front of the barrier. 

Unlike the case of a wide barrier, in the case of the 
8 function nonlinear barrier potential, localized steady 
states exist only at v < v c < v s where v s is the sound 
velocity. In our case v s = 1. Critical velocity v c is deter- 
mined by the potential strength 7 



7 



16(1 - v 2 ) 2 (2v 2 - 3 + ajvc)) 1 / 2 
(6v 2 - 3 + a(v c )) 2 {-2v 2 - 1 + a{v c )y/ 2 ' 



(24) 



where a(v c ) = \/9 — ^v 2 + 4w^. 

In order to cover a wide range of velocities we have 
carried out numerical simulations of the flow of a BEC 
through the delta potential nonlinear barrier moving with 
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FIG. 6: Time evolution of steady solutions with the param- 
eter xq = xo+ = 0.752048 and xo = xo- = 0.350966 corre- 
sponding relatively to (a) stable and (b) unstable BEC flows 
past a nonlinear repulsive delta potential barrier in subcriti- 
cal regime. The other parameters v — 0.65 and 7 = 0.5. For 
this case v c = 0.663946. 

small acceleration beginning from zero velocity. Fig.[7]dc- 
picts the time evolution of a BEC flow when the acceler- 
ation a = 0.004. The barrier potential strength 7 = 0.5. 
The initial wave packet is taken in the form of Eq. (|2ip . 

Time interval 0<f<165(0<v< v cr ) corresponds to 
a superfluid flow. At tmes 170 < t < 250 (v cr < v < v s ) 
one can observe generation of grey solitons chain. In 
time interval 250 < t (v s < v) corresponding to trans- 
critical flow of a BEC at supersonic velocity one can ob- 
serve qualitatively the same wave pattern obtained in the 
work 1] where a dispersive shock propagates upstream 
with generation of soliton-like waves propagating down- 
stream. 



III. CONCLUSION 

In conclusion, we studied steady flow in a defocusing 
quasi ID BEC moving through a nonlinear repulsive bar- 
rier. Such a kind of barriers can be formed by variation 
of the atomic scattering length of BEC in space. For the 
case of a wide nonlinear barrier we have found critical 



FIG. 7: Time evolution of a BEC flow through the delta po- 
tential nonlinear barrier moving with the acceleration a = 
0.004. The barrier potential strength 7 = 0.5, initial velocity 
of the flow vo — 0. The initial wave packet is taken in the 
form Eq. (pTj) . 



velocities of steady flows. Within the interval of veloci- 
ties v_ < v < v + , in the transcritical regime we observed 
generation of a slow moving train of dark solitons. At ve- 
locities above supercritical the train disappears. At the 
same time in this regime one can observe formation of a 
hump localized at the place of the barrier. 

For the case of a 6 function nonlinear barrier potential 
the dependence of the steady solution parameters and a 
critical velocity on the potential strength 7 was found in 
analytical form. As numerical simulations show, in sub- 
critical regime v < v c an unstable solution decays into a 
gray soliton moving upstream and a stable solution local- 
ized at the barrier position. The decay is accompanied by 
a dispersive shock wavepropagating downstream in front 
of the barrier. 

The dynamics of flows past through a linear and non- 
linear barriers are qualitatively similar except the follow- 
ing. In the case of a wide linear barrier, the superfluidity 
is broken at any small velocities if the barrier potential 
strength greater than some threshold value (see Fig. 2 in 
1]). For a wide nonlinear barrier an interval of velocities 
< v < u_ always exists, where the flow is superfluid 
regardless of the barrier potential strength. 

When using the optically induced Feshbach resonance 
technique to generate a repulsive nonlinear barrier by 
focused laser beam, one should in general take into ac- 
count the losses, induced by spontaneous emission of 
atoms. Phenomenologically it can be described by adding 
a nonlinear loss term —ij\u\ 2 u in the GP equation. Atom 
feeding can be described by linear gain term ictu. This 
case requires a separate investigation. It should be noted 
that this problem relates to one considered in the recent 
work [l5j], where the flow of polariton condensate [l6| 
past a linear barrier was studied taking into account lin- 
ear amplification and nonlinear damping. 



6 



Acknowledgments No.PIIF-GA-2009-236099(NOMATOS). 

Authors are grateful to E.N. Tsoy for fruitful dis- 
cussions. F.Kh.A. acknowledges a Marie Curie Grant 



[1] A.M. Leszczyszyn, G.A. El, Yu.G. Gladush, and A.M. 

Kamchatnov, Phys.Rev. A 79, 063608 (2009). 
[2] GE. Astrakharchik and L.P. Pitaevskii, Phys.Rev. A 70, 

013608 (2004). 

[3] P.Engels and C. Atherton, Phys.Rev.Lett. 99,160405 
(2007). 

[4] B. Malomed and M. Azbel, Phys.Rev. B 47, 10402 
(1993). 

[5] V. Hakim, Phys.Rev. E 55, 2835 (1997). 
[6] L. P. Pitaevskii and S. Stringari, Phys. Lett. A 235, 398 
(1997). 

[7] R.P. Feynman, in Progress in Low Temperature Physics, 
edited by C.J. Coster (North-Holland, Amsterdam, 
1955), Vol.1, p.17. 

[8] C. K. Law, C. M. Chan, P. T. Leung, and M.-C. Chu, 



Phys. Rev. Lett. 85, 1598 (2000). 
[9] M. Haddad and V. Hakim, Phys. Rev. Lett 87, 218901 
(2001). 

[10] N. Pavloff, Phys.Rev. A 66, 013610 (2002). 

[11] S. Inouye et ah, Nature 392, 151 (1998). 

[12] P.O. Fedichev, Yu. Kagan, GV. Schlyapnikov, and 

J.T.M. Walraven, Phys.Rev.Lett. 77, 2913 (1996). 
[13] R. Yamazaki, S. Taie, S. Sugawa, and Y. Takahashi, 

Phys. Rev. Lett. 105, 050405 (2010). 

[14] N. Berloff and V.M. Perez-Garcia, I arXiv: 1006.44261 
[15] A.M. Kamchatnov and Y.V. Kartashov, EPL 97, 1006 

(2012). 

[16] A. Amo et al, Science, 332, 1167 (2011). 



